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Abstract 

This article focuses on techniques for position reconstruction of acoustic point sources with the AMADEUS setup 
consisting of 36 acoustic sensors in the Mediterranean Sea. The direction reconstruction of an acoustic point source 
utilizes the information of the 6 small-volume hydrophone clusters of AMADEUS individually. Source position 
reconstruction is then done by combining the directional information of each cluster. The algorithms for direction 
and position reconstruction are explained and demonstrated using data taken in the deep sea. 
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1. Introduction 

AMADEUS CI (Antares Modules for Acoustic 
Detection Under the Sea) is a system for acoustic 
detection within the ANTARES [2] (Astronomy 
with a Neutrino Telescope and Abyss environ- 
mental RESsearch) infrastructure. The former is 
designed to investigate the feasibility of acous- 
tic detection of ultra high energy neutrinos. The 
ANTARES experiment is located 25 km off the 
southern French coast at a water depth of 2500 m 
and consists of 12 vertical support structures, so 
called lines, of 480 m length. They are anchored to 
the sea floor over an area of about (180 x 180) m 2 , 
held vertically by buoys and support a total of 
~ 900 photo-multipliers. An additional line, the 
Instrumentation Line, is used for environmental 
monitoring and multidisciplinary research. The 
AMADEUS system consists of six local clusters of 
six sensors each. Three of these clusters are located 
on the 12th line while the other three clusters are 
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located on the Instrumentation Line (IL07), see Fig. 
[TJ AMADEUS was designed to investigate short 



Acoustic Storey 
( Pointing Down ) 



Acoustic Storey 
(Standard) 




Cable to shore 



Fig. 1. The ANTARES detector. The AMADEUS clusters are 
located on the Instrumentation Line (IL07) and Line 12 (L12). 

and long distance correlations, study the deep 
sea noise, signal filtering and source positions re- 
construction techniques, all necessary to evaluate 
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the feasibility of acoustic particle detection with 
a possible future large volume detector. Using lo- 
cal clusters allows for correlation studies ranging 
from about 1 m to 350 m, and also offers the ad- 
vantage of quick direction reconstruction, thereby 
simplifying source position reconstruction. 

2. Direction reconstruction 

Each of the six local clusters of AMADEUS com- 
prises six hydrophones. The hydrophones within a 
local cluster have typical distances of 1 m between 
the sensors and thus a transit time of < 1 ms for 
an acoustic wave front. This allows for quick di- 
rection reconstruction of a passing pressure wave. 
Two methods for direction reconstruction with the 
individual clusters will be discussed. One of the 
methods is beamforming and the other one uses 
the time differences of the detected signals in the 
individual acoustic sensors of a storey. 

2.1. Beamforming 

This method requires the knowledge of the hy- 
drophone positions, within a cluster, as well as 
synchronized data. Beamforming is done by creat- 
ing a sound intensity plot scanning all directions 
in space (4zr) for the actual source direction. Given 
the six hydrophone coordinates f n (n = 1,2..., 6), 
the signal p n of every hydrophone will be shifted 
in time corresponding to the difference in path 
length of the sound wave reaching the respective 
hydrophone. Hence every direction in space corre- 
sponds to a set of time differences At n in the data. 
For a direction fc, the beam-forming output at time 
t is given by 



N 



(1) 



n=l 



where w n represent weighting factors for every in- 
dividual hydrophone. These factors are adjusted 
to match the directional sensitivity of each hy- 
drophones. For the following calculations w n = 1 
was used. The time differences At n are computed 
assuming a plane wav^H The algorithm scans 4zr 
with a predefined step size in Q, (p by applying the 
calculated time differences to the data, assuming 



1 The assumption of a plane wave is a sufficiently good ap- 
proximation if the distances are large compared to the dimen- 
sions of the hydrophone antenna. 



a constant speed of sound. Once all directions are 
scanned, the maximum value of the produced out- 
put identifies the direction of the signal. For this 
direction all signals are shifted such that they add 
up constructively. The application of the algorithm 
increases the signal to noise ratio for the direc- 
tion the signal comes from, creating the mentioned 
maximum. Figure |3] shows the output of the beam- 
forming algorithm applied to the signal in Fig. 12 
The beam-forming output (Fig. [3j> shows a well 
defined maximum at = 59° and cp = —29°. The 
error is governed by the binning of the algorithm 
and not by the algorithm itself. Local maxima and 
the visible patterns in Fig. [3] correspond to direc- 
tions where the time shift results in the construc- 
tive addition of two or more signals. 
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Fig. 2. Signals recorded by one local cluster, consisting of six 
hydrophones. 
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Fig. 3. Output of the beam-forming algorithm applied to the 
signal shown in Fig. |2 The maximum, indicating the signal 
direction, can be identified at 6 = 59° and (p — —29°. 
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The advantage of this method is that it can be 
applied to untriggered data. A disadvantage of this 
method is that it is rather time consuming. 

2.2. Time difference method 

Another method for directional reconstruction 
is using the time differences between signals, de- 
tected by the hydrophones in a cluster, of a passing 
pressure wave. This method requires, in addition 
to the knowledge of the hydrophone position and 
synchronized data, a signal identification. Once a 
signal is identified, e.g. by a simple threshold trig- 
ger for large signal amplitudes, the signal arrival 
times are used for the direction reconstruction. Di- 
rection reconstruction is done by comparing the 
measured signal times to expected signal arrival 
times. The expected time t nexpected (6, q>) for an ar- 
riving wave front is taken from a lookup table and 
subtracted from the measured time. That is done 
for 4 n in the desired angular resolution, where the 
minimum indicates the reconstructed direction: 



min 



6 / 



t„, 
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expected 



(2) 



This technique is similar to beam-forming, but 
much faster. Applied to the data in Fig. |2] the algo- 
rithm produces the r 2 data shown in Fig. H] and 
A clear minimum can be identified in both and (p. 
The values match the ones obtained by the beam- 
forming algorithm. The error is again governed by 
the binning of the algorithm. 





Fig. 5. Cut through the 6 — cp plane at <p m i n of the minimum 
computed with the time difference method. 



Though the time difference method is much 
faster, beam-forming has the advantage of not 
having to identify a signal. 

Figure shows a qualitative mapping of the ar- 
rival directions of transient acoustic signals origi- 
nating in the surrounding of the ANTARES detec- 
tor. The direction reconstruction was done using 
the time difference method. The data is from the 
time period 29.01.2008 to 05.06.200© using a min- 
imum bias trigger which is taking unfiltered data 
for 10 seconds every 30 minutes. Taking this data 
set, the directions of all included signals with an 
amplitude greater than 8 -RMS (about 200000 sig- 
nals) were reconstructed. The RMS was calculated 
for a 0.1 s time interval around the signal time. 
A virtual observer residing on the second cluster 
from the bottom of the IL07, about 180 m above the 
sea bed, looking westward towards the horizon of 
the storey, is defining the origin of the coordinate 
system, see Fig. Thus 90° (-90°) in longitude (<p) 
is north (south), while 90° (-90°) in latitude (0) is 
looking upwards (downwards). As expected, the 
majority of the signals is coming from above the 
detector, including all kinds of transient signals, 
as e.g. produced by dolphins. In the lower hemi- 
sphere a few sources can be seen amongst which 
a clear structure can be identified. The Signals be- 
tween 0° and -120° in <p and -30° to -70° in 
originate from the ANTARES acoustic position- 
ing system, emitting signals at the bottom of each 
string. 



Fig. 4. Cut through the 6 — (p plane at 9 m { n of the minimum 
computed with the time difference method. 



2 Most of the data was taken when the detector consisted of 
10 lines. 
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Fig. 6. A qualitative mapping of the arrival directions of tran- 
sient acoustic signals. The majority of the signals come from 
above. In the lower hemisphere the ANTARES acoustic posi- 
tion system can be identified. 

3. Source reconstruction 



For a given local cluster such as in the 
AMADEUS setup, directions are the prime ob- 
servables of acoustic signals. The identified di- 
rections then allow for a simple source position 
reconstruction. After the successful direction re- 
construction, each antenna i that detected a signal 
will point to a direction k\ identifying where an 
event came from. Since the cluster positions S\ are 
known, one obtains a line pointing to the event 
for each cluster, 



di = Si + rijkj , rii e IR. 



(3) 



The source position of a signal is theoretically 
given by the intersection point of the lines dj. Due 
to inherent uncertainties in the direction recon- 
struction, mainly due to the one degree grid, the 
source location is where the lines are closest to 
each other. Localising the point of closest approach 
for all lines is realized by calculating the square 
distance from a point s to the reconstructed lines, 

L?(s) = (V- (Sf + ft;fc/)) witn n i = (s - Si) • ki , (4) 

and minimising the sum of these distances. 

Since the cluster positioning and orientation was 
not yet implemented, a full source position recon- 
struction was not possible at the time of this anal- 
ysis. However, as the tilts of the local clusters, in- 
duced by the sea current, are small, in the order of 
one degree, it is possible to reconstruct the 6 an- 
gle with that uncertainty. The circles in Fig. rep- 
resent the possible solutions of a source given by 



the reconstructed angle and the distance of each 
storey to the sea floor. The signal origin can clearly 
be identified as the acoustic positioning emitter of 
the 7th line (L7), as can also be deduced from the 
signal properties. 
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Fig. 7. Layout of ANTARES on the sea floor in UTM coordi- 
nates. Each point indicates the position of a line. The positions 
of the Instrumentation Line (IL07), Line 12 (L12) and Line 7 
(L7) are labeled. 



4. Summary 

Techniques for source reconstruction using lo- 
calized clusters have been presented. Source lo- 
cation is done by using the local clusters for di- 
rection reconstruction. For direction reconstruction 
two methods were investigated. One is a beam- 
forming algorithm only requiring the hydrophone 
positions, within a local cluster, and synchronized 
data. The other method additionally requires a sig- 
nal identification and uses time differences. Al- 
though both methods work well, using time dif- 
ferences proved to be much faster. A technique for 
source reconstruction, based on the directional in- 
formation of the local clusters, was explained and 
its potential demonstrated. 

This study was supported by the German govern- 
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